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Abstract 

The universal primer three-primer approach can dramatically reduce the cost when genotyping the microsatellites. One 
former research reported four universal primers that can be used in singleplex and multiplex genotyping. In this study, we 
proposed an alternative suite of universal primers with four dyes for genotyping 8-12 loci in one single run. This multiplex 
method was tested on Tetranychus truncatus. Published microsatellite loci of T. kanzawai, Frankliniella occidentalis and 
Nilaparvata lugens were modified as needed and also tested. The robustness of the method was confirmed by comparing 
with singleplex using multiple fluorophores and genotyping two populations of T. truncatus. This method showed lower 
signal strength than the singleplex three-primer system, but it was still sufficient to determine the fragment length. The cost 
of such a project can be reduced dramatically when many loci of different species are involved. In this way, laboratories 
performing population genetic analyses or studying several different species may benefit from the use of this cost-effective 
protocol. 
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Introduction 

Microsatellites were first genotyped by Litt and Luty [1]. There 
are now a number of methods available for the identification of 
SSR alleles. Methods based on radioactive isotopes and silver 
staining have been almost completely replaced by the more rapid 
and precise capillary electrophoresis fluorescence detection system, 
which has revolutionized the field. Some 90 out of 100 articles 
relying on SSRs published in 2009 and 2010 in the journal 
Molecular Ecology involved subjects genotyped using this automatic 
capillary electrophoresis system [2]. 

In capillary electrophoresis systems, forward or reverse primer of 
one specific locus was usually directly labeled with a fluorescent tag. 
Then this fluorescence is incorporated into PGR products, which 
can be detected in the capillary sequencer. Then the length of the 
PGR products can be analyzed. However, the use of this method is 
limited because the fluorescent dyes are very expensive and the cost 
increases considerably when a large number of loci are involved. 
Meanwhile, an economical approach was proposed and improved 
[3-7], which have been cited for more than 1300 times. The 
universal primer three-primer approach is based on a universal 
fluorescent primer (M 1 3) which is identical to the sequence added at 
5 ' end of each forward primer. During the latter PGR cycles, this 
universal fluorescent primer hybridizes with the complementary 
region generated in the early PGR cycles. In this way, instead of 
synthesizing fluorescent-labeled primer for each SSR locus, only 
one universal labeled primer must be fluorescent. 

Theoretically, this dye-labeled universal primer method can be 
multiplexed [6]. However, the universal primer three-primer 



approach limits multiplexing if loci sizes do not overlap. Although 
the Ml 3 sequence is the most widely used tail, any sequence that is 
not homologous to the target genome could be used. Missiaggia 
and Grattapaglia tested Eucalyputs DNA samples using three 
universal tags DYS437, D8S1132 and D12S1090, which are 
derived from human microsatellites [8]. Moreover, Blacket et al. 
demonstrated that four universal primers can be combined with 
four fluorophores to co-amplify multiple loci via multiplex PGR 
[9]. This strategy, called strategy Ill-Multiplex with multiple 
fluorophores, is highly cost-effective method and has been proven 
to be feasible [10-13]. 

Here the concept and methodology of strategy III described by 
Blacket et al were kept [9] . The simultaneous use of another four 
universal primers combined with four fluorophores to genotype 8- 
12 loci in a single PGR run is proposed here. The effectiveness of 
the new set of universal primers in multiplex genotyping was tested 
for the detection of microsatellites of spider mites Tetranychus 
truncatus and further validated by genotyping two populations of T. 
truncatus. Moreover, published microsatellite loci of T. kanzawai, 
Frankliniella occidentalis and Nilaparvata lugens were modified as 
needed for universality test. 

Materials and Methods 

Mite collection and identification 

Mites were collected in 2011 and 2012 from five provinces in 
Ghina (Table 1). The mites were reared on detached bean leaves 
[Phaseolus vulgaris) at 25°G in a 16 h light/8 h dark cycle. 
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Table 1. Collection Information for samples used In this study. 



Population name 


Site 


Location 


Species 


Sampling date 


No. individuals used 


Host 


BC 


Baotou 


40°39'N,109°50'E 


T. truncatus 


Sep. 2011 


4 


Phaseolus vulgaris 


SQ 


Shenyang 


4r49'N, 123°33'E 


T. truncatus 


Sep. 2011 


50 


Solanum melongena 


XZ 


Xuzhou 


34°44'N,116°56'E 


T. truncatus 


Jul. 2012 


42 


Gossypium hirsutum 


XN 


Xining 


36°29'N,101°34'E 


T. kanzawai 


Sep. 2011 


4 


Livistona cliinensis 


PZ 


Pengze 


29°51'N,116°35'E 


T. kanzawai 


Jul. 2012 


4 


G. hirsutum 


HB 


Harbin 


45°45'N,126°38'E 


F. occidentalis 


Jul. 2010 


4 


Tagetes erecta 


DL 


Dali 


25°36'N,100°15'E 


F. occidentalis 


Jul. 2011 


4 


Canna indica 


SY 


Sanya 


29°5rN,116 35'E 


N. lugens 


Aug. 2013 


4 


Oryza sativa 


GL 


Guilin 


29°51'N,116°35'E 


N. lugens 


Aug. 2013 


4 


0. sativa 



doi:10.1371/journal.pone.0092826.t001 



The identification of spider mites has long been an issue because 
different species are very similar. In one previous study, both a 
morphological approach based on the shape of the aedeagus and a 
molecular approach based on the internal transcribed spacer (ITS) 
region sequence were used to confirm the identities of spider mite 
individuals, which proved to be accurate and efficient [14]. This 
method was also used in the present study. 

No specific permits were required for the study, (a) No specific 
permissions were required for the mite collections because the 
mites are pests on crops; (b) The location is not privately-owned; 
(c) The study did not involve endangered or protected species. 

DNA extraction and primer preparation 

Female T. kanzawai were direcdy picked from detached bean 
leaves and then washed with MiUi-Q^ water. Genomic DNA was 



extracted using a DNeasy tissue kit (Qiagen, U.S.) following the 
protocol described by Tsagkarakou et al. [15]. The DNA of T. 
truncates was extracted in 20 1 2 using the same protocol and stored 
in — 20°C. The DNA samples of occidentalis and M. lugens were 
provided by our colleagues (collection information see Table 1). 
The quality of all DNA samples was estimated using a NanoDrop 
ND-1000 spectrophotometer. 

Twelve primer parrs (Table 2) were selected for amplification of 
the SSR loci of T. truncates with the help of Multiplex Manager 1 .2 
[14,16]. Each specific forward primer was connected with one of 
the 5' universal primer sequence tails which were T7 (5'- 
TAATACGACTCACTATAGGG), Ml 3 (5'-TGTAAAAC- 
GACGGCCAGT), M13R (,5'-CAGGAAACAGCTATGACC) 
or AP2 (5'-CTATAGGGCACGCGTGGT) (Table 2). The stand 
alone universal primer Ml 3, AP2, M13R and T7 were 



Table 2. Universal primer, dye and group of microsatellite primers of Tetranychus species used in tliis study. 



Species Locus(Ps) Universal primer added to forward primer Dye Group 



T. truncatus 


TUFG87 


M13 


FAM 


Pa 




TUFG57 


AP2 


VIC 






TUFG127 


M13R 


NED 






TUFG130N 


T7 


PET 




T. truncatus 


TUFG177 


M13 


FAM 


Pb 




TUFG55 


AP2 


VIC 






TUFG146 


M13R 


NED 






TUFG147N 


T7 


PET 




T. truncatus 


TUFG33 


M13 


FAM 


Pc 




TUFG43 


AP2 


VIC 






TUFG67 


M13R 


NED 






TUFG158 


T7 


PET 




T. kanzawai 


TklVlS002 


Ml 3 


FAM 


Pe 




TkMSOlO 


AP2 


VIC 






TkMSOl 1 


M13R 


NED 






TkMS014 


T7 


PET 




T. kanzawai 


TkMSOOe 


Ml 3 


FAM 


Pf 




TkMSOl 3 


AP2 


VIC 






TkMSOl 5 


M13R 


NED 





Pa+b/ all the PjS in primer mix Pa and Pb were pooled; Pa+b+c- all the P^s in primer mix P^, Pb and Pc were pooled. 
doi:10.1371/journal.pone.0092826.t002 



PLOS ONE I www.plosone.org 



2 



March 2014 | Volume 9 | Issue 3 | e92826 



Genotyping by Universal Primers in IVlultiplex 
Table 3. Universal primer and dye of microsatellite primers of F. occidentalis and N. lugens used in this study. 



Species 


Locus 


Universal primer* Dye 


Species 


Locus 


Universal primer* 


Dye 


F. occidentalis 


WFT20 


M13 


FAM 


N. lugens 


NL22 


M13R 


NED 




WFT28 


T7 


PET 




NL41 


M13 


FAM 




WFT37 


MIS 


FAM 




NL74 


M13R 


NED 




WFT51 


Ml 3 


FAM 




NL121 


T7 


PET 




WFT64 


AP2 


VIC 




NL122 


AP2 


VIC 




WFT66 


M13R 


NED 




NL140 


M13 


FAM 




WFT83 


T7 


PET 




NL157 


AP2 


VIC 




WFT98 


T7 


PET 




NL158 


AP2 


VIC 




WFT104 


M13R 


NED 




NL161 


T7 


PET 




WFT108 


AP2 


VIC 




NL162 


T7 


PET 












NL167 


M13 


FAM 












NL177 


M13R 


NED 



*Universal primer added to corresponding forward primer. 
doi:l 0.1 371 /journal.pone.0092826.t003 



fluorescent-labeled with FAM, VIC, NED and PET respectively. 
All primers were ordered from Life Technologies and dissolved in 
TE (10 mM Tris, 1 mM EDTA, pH 8.0) to produce a 20 ^iM 
stock solution for fluorescent-labeled universal primer and locus- 
specific reverse primer and a 5 |lM stock solution of locus-specific 
forward primer with universal tail (moles of labeled forward 
primer: reverse primer: dye-labeled universal primer = 1:4:4 [6]). 



Primer mix is here designated Ps (three primers were mixed for 
each specific locus and are collectively called Ps, as in PxuFGsy)- 
Pa, Pb, Pg, Pa+b and Pa+b+c were also prepared. They contained 
all corresponding P^s (see Table 2) at equal volume and then stored 
in small aliquots at — 20°C in the dark. 




160 180 200 220 240 260 280 300 320 
1200 i-H 1 1 1 1 1 ' : 1 : ' 1 \ K- 



900- 



SOO- 




Figure 1. Multiplex PCR results of T. truncatus BC02. Three panels refer to singleplex, samples amplified by and Pb separately then were 
pooled for capillary loading and samples amplified by Pa+b respectively. PI: TUFG130N. VI: TUFG55. F1: TUFG177. N1: TUFG127. V2: TUFG57. N2: 
TUFG146. P2: TUFG147N. F2: TUFG87. 
doi:1 0.1 371 /journal.pone.0092826.g001 



PLOS ONE I www.plosone.org 



3 



March 2014 | Volume 9 | Issue 3 | e92826 



Genotyping by Universal Primers in IVlultiplex 




Figure 2. Multiplex PCR results of F. occidentalis HB04. Two panels refer to singleplex and sample amplified by ten loci in one multiplex PCR 
respectively. PI: WFT28. N1: WFT66. VI: WFT108. N2: WFT104. F1: WFT20. P2: WFT98. V2: WFT83. F2: WFT83. V3: WFT64. F3: WFT51. 
doi:1 0.1 371 /journal.pone.0092826.g002 



Genotyping 

The usefulness of these four universal primers was tested with 
eight T. trmcatus samples from population XZ and BC in three 
levels: singleplex PCR, multiplex PCR in which each dye labeled 
one locus and the multiplex PCR in which each dye labeled more 
than one locus. At level I, PCR was carried out in 5 |iL volumes 
using an Applied Biosystems Veriti Thermal Cycler (Applied 
Biosystems). Each assay contained 2.1 |iL RNase-free water, 1 |iL 
template DNA, 0.15 |iL Pg and 1.75 |xL Type-it Microsatellite 
PCR Master Mix (2X, Qiagen). At level II and III, PCR was 
performed in a final volume of 5 |a,L containing 1 .95 |.iL RNase- 
free water, 1 template DNA, 0.3 nL Pa/Pb (level II) or Pa+b/ 
Pa+b+c (level III) and 1.75 nL Type-it Microsatellite PCR Master 
Mix (2X, Qiagen). At all three levels, PCR was performed on the 
same DNA samples in the same 96 well plate simultaneously under 
the same cycling conditions: an initial activation step for 5 min at 
95°C; 8 cycles of 95°C for 30 s, 57°C for 90 s, 72°C for 30 s, 
followed by 14 cycles of 95°C for 30 s, annealing at 57-50°C 
decreased by 0.5°C per cycle for 90 s, 72°C for 30 s, and then 12 
cycles of 95°C for 30 s, 52°C for 90 s, 72°C for 30 s and a final 
extension at 60°C for 30 min. 

At level I, PCR products of the same DNA samples were all 
mixed and pooled within a well. The 1 .0 |iL of the mixture of level 
I or 1.0 nL PCR product in level II or III was added to 10 |iL 
formamide and then electrophoresed on the ABI 3130 capillary 
sequencer (Applied Biosystems) along with the GeneScan-500 LIZ 
size standard. The data were collected automatically and alleles 
were called using GeneMapper v4.0 (Applied Biosystems). 



Method further validation 

The usefulness of multiplex PCR with these four alternative 
universal primers was validated through the genotyping of two 
populations of T. trmcatus and microsateUites of other species 
which have been pubUshed. Two populations of T. truncatus were 
amplified by primer mix Pa+b- The PCR procedure and condition 
were the same as above. For further validation, seven microsat- 
eUites primers that used in T. kanzawai were chosen [17]. In the 
original paper, one primer in each locus was modified with 
fluorescent dye 6-FAM, HEX, or TAMRA. Therefore, these 
primers were directly modified as shown in Table 2 to 
accommodate our protocol. PCR mixtures were prepared as 
above with the following profde: 95°C for 5 min; 35 cycles of 95°C 
for 30 s, 52°C for 90 s, and 72°C for 30 s; and a final elongation at 
60°C for 30 min. The genotyping procedure was the same as of 
the one used for T. truncates. 

Further, ten and twelve microsatellite loci with null allele 
frequency<0.1 were selected for genotyping eight samples of 
occidentalis [18] and jV. lugens [19] respectively. The primers of each 
locus were modified as shown in Table 3. PCR mixtures were 
prepared and genotyped as above. PCR profile for F. occidentalis 
was 5 min at 95°C, 32 cycles of 95°C for 30 s, 52°C for 90 s, 72°C 
for 30 s, followed by a final extension at 60°C for 30 min. PCR 
profile for jV. lugens was 5 min at 95°C, 12 cycles of 95°C for 30 s, 
annealing at 58-52°C decreased by 0.5°C per cycle for 90 s, 72°C 
for 30 s, followed by 20 cycles of 95°C for 30 s,52°C for 90 s, 
72°C for 30 s, and then 60°C for 30 min. 
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Figure 4. Microsatellite profiles of two DNA samples amplified by loci Ti<MS010 (left peai<s) and TkMS013 (right peaks) through 
three methods generated on 3130 capillary sequencer, a) DNA sample PZ02; b) DNA sample PZ04. Each part of the figure contains three 
screenshots which referred to singleplex, samples amplified by Pe and Pp then pooled for capillary loading and samples amplified by Pe+f using the 
method described herein. 
doi:1 0.1 371/journal.pone.0092826.g004 



Data analysis 

Only genotyping data from populations of T. truncak.s were 
analyzed as follows. NuU allele frequencies were determined with 
Micro-Checker version 2.2.3 with the Brookfield's nuU aUele 
estimator 2 [20,21]. Genotype errors such as stuttering or large 
allelic dropout were also detected using Micro-Checker. Genepop 
4.0.10 was used to search for deviations from Hardy-Weinberg 
equUibrium (HWE) at each locus/poptilation combination using 
Fisher's exact tests [22]. The population genetic diversity 
parameters of the population such as total alleles per locus, 
observed heterozygosity and expected heterozygosity were calcu- 
lated in GenAlEx 6.5 [23]. Polymorphism information content 
(PIC) was calculated using Cervus version 3.0 and inbreeding 
coeflTicient (i?is) with the FSTAT 2.9.3.2 [24,25]. 

Results and Discussion 

Here the methodology described by Schuelke and Blacket et al. 
were inherited and extended by proposing another suite of 
universal primers for simultaneously genotyping 8-12 loci in a 
single run. [6,9]. In addition, direct comparisons between 
singleplex and multiplex PGR were presented and different levels 
of multiplex genotyping were evaluated. This significantly reduces 
the time and cost required for medium-throughput microsatellite 
genotyping experiments. Twelve SSR markers from a previous 



paper served as volunteers to genotype T. truncatus DNA samples. 
The universality of the method was tested by modifying existing 
primers of microsatellite loci then genotyping T. kanzawai, F. 
occidmtalis and jV. lugens. The method was tested further by 
genotyping two populations of T. truncatus. 

PGR products of the same DNA samples were grouped in one 
project to evaluate the quality of allele peaks generated using this 
method (two multiplex results see Fig 1 and Fig 2, others results of 
Tetmnjchus species can be found in Figures S1-S38 in Online 
Resources SI, results of lugens and F. occidentalis can be found in 
Figures S39-S42 and Figures S43-S46, respectively in Online 
Resources S2). As expected, the majority of the signals were 
weaker than the conventional singleplex PGR because all pairs of 
primers competed for a certain amount of resources, such as 
templates and enzyme. However, the signals strength was still 
sufficient to measure the fragment length. The expected micro- 
satellite morphology was observed. No difference in fragment 
length was observed. For the T. truncatus, all the loci were amplified 
well and showed the same allele size, which was sufficient for 
studies of population genetics. The DNA sample XZ04 was 
amplified with loci TUFG55, TUFG57, TUFG87 and TUFG177 
using all three methods described above. Results are shown 
together in Fig. 3 (other results including eight samples of XZ 
amplified by Pa+b+c available in Figures S25 and S26 in 
Online Resources SI). For T. kanzawai, locus TkMS013 was 
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unable to amplify any DNA samples using any of the normal 
methods of singleplex PGR. It did work, however, when the 
multiplex method was used (Fig 4). This result persisted in several 
trials. This success cannot be attributed to the primers because the 
primer mixes for multiplex came from primer mix Ps (three- 
primer system for each specific locus) for a single PGR. All the 
other loci of T. kanzawai performed well for all three methods. 
Another phenomenon was observed at some loci as shown in locus 
TkMSOlO (Fig 4). These loci showed greater signal of stutter peak 
which adjacent to the true peak. This variation was stable 
throughout the experiments. Thus, once bins were determined (no 
matter the stutter peak or the true peak), the alleles calling 
remained consistent. 

Examination of genotyping quality by two populations of T. 
tmncatus showed the results to be consistent with those of a previous 
paper on SSR development [14]. As expected, 2 bp differences in 
allele size were observed in the loci TUFG130N and TUFG147N 
as T7 tag was used instead of Ml 3. The quality of the genotype 
data was confirmed using Micro-Ghecker 2.2.3. It did not show 
stuttering or large allelic dropout and the characteristics of loci 
were consistent with the results of a previous paper and of a 
structural study that is currendy in progress (data not shown). For 
these eight loci in Pa+b, the PIG ranged from 0.402 (TUFG146) to 
0.790 (TUFG87). The mean number of alleles per locus was 6.5. 
The observed and expected heterozygosity across the two 
populations ranged from 0.048 (TUFG146 in XZ) to 0.820 
(TUFG57 in SQ) and from 0.091 (TUFG146 in XZ) to 0.752 
(TUFG57 in SQ), respectively. Loci TUFG57, TUFG87 and 
TUFG177 in population XZ deviated from HWE due to 
heterozygote deficits with null allele frequency >0.120 (0.120- 
0.230) and inbreeding coeflTicient >0.206 (0.206-0.364). No 
departure from HWE was observed at any of locus in the SQ, 
population and there was no indication of any null allele. AU 
inbreeding coefficients were <0.133. AU these results demonstrate 
that all the amplifications took place in a correct and consistent 
manner when the multiplex method was used. 

This multiplex protocol reduces the cost dramatically because it 
makes full use of the fluorescent dyes. Gosts can also be decreased 
by using Qiagen Multiplex PGR Kits sparingly. Here the PGR 
final volume was reduced to 5 (iL with a buffer concentration of 
0.7 X, without compromising performance or specificity [26]. The 
PGR products were a mixture of several amplifications, ^vliicli 
could be diluted and directiy mounted on the Sequencer. Here, we 
presented another suite of universal primers with lower annealing 
temperature (approximately 52°G) facilitating the modification of 
existing SSR primers which have annealing temperatures between 
52°G and 58°C. By directly comparing the different level of 
multiplex, we found that most loci in multiplex showed lower 
signal than in singleplex, which can be reasoned by the 
competition between loci. However, no significant signal reduction 
was found between level II and level III. At level II, one locus 
monopolized one fluorescent dye while at level III one fluorescent 
dye were shared by at least two loci. Thus multiplex with multiple 
fluorescent dyes is preferred. 

For optimized results the following must be held true: (1) the 
amount of labeled forward primer: reverse primer: dye-labeled 
universal primer = 1:4:4. This ratio must hold true for each 
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